home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / sggbak.z / sggbak
Text File  |  1998-10-30  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGGGGGBBBBAAAAKKKK((((3333FFFF))))                                                          SSSSGGGGGGGGBBBBAAAAKKKK((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGGBAK - form the right or left eigenvectors of a real generalized
  10.      eigenvalue problem A*x = lambda*B*x, by backward transformation on the
  11.      computed eigenvectors of the balanced pair of matrices output by SGGBAL
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SGGBAK( JOB, SIDE, N, ILO, IHI, LSCALE, RSCALE, M, V, LDV,
  15.                         INFO )
  16.  
  17.          CHARACTER      JOB, SIDE
  18.  
  19.          INTEGER        IHI, ILO, INFO, LDV, M, N
  20.  
  21.          REAL           LSCALE( * ), RSCALE( * ), V( LDV, * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      SGGBAK forms the right or left eigenvectors of a real generalized
  25.      eigenvalue problem A*x = lambda*B*x, by backward transformation on the
  26.      computed eigenvectors of the balanced pair of matrices output by SGGBAL.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      JOB     (input) CHARACTER*1
  31.              Specifies the type of backward transformation required:
  32.              = 'N':  do nothing, return immediately;
  33.              = 'P':  do backward transformation for permutation only;
  34.              = 'S':  do backward transformation for scaling only;
  35.              = 'B':  do backward transformations for both permutation and
  36.              scaling.  JOB must be the same as the argument JOB supplied to
  37.              SGGBAL.
  38.  
  39.      SIDE    (input) CHARACTER*1
  40.              = 'R':  V contains right eigenvectors;
  41.              = 'L':  V contains left eigenvectors.
  42.  
  43.      N       (input) INTEGER
  44.              The number of rows of the matrix V.  N >= 0.
  45.  
  46.      ILO     (input) INTEGER
  47.              IHI     (input) INTEGER The integers ILO and IHI determined by
  48.              SGGBAL.  1 <= ILO <= IHI <= N, if N > 0; ILO=1 and IHI=0, if N=0.
  49.  
  50.      LSCALE  (input) REAL array, dimension (N)
  51.              Details of the permutations and/or scaling factors applied to the
  52.              left side of A and B, as returned by SGGBAL.
  53.  
  54.      RSCALE  (input) REAL array, dimension (N)
  55.              Details of the permutations and/or scaling factors applied to the
  56.              right side of A and B, as returned by SGGBAL.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSGGGGGGGGBBBBAAAAKKKK((((3333FFFF))))                                                          SSSSGGGGGGGGBBBBAAAAKKKK((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      M       (input) INTEGER
  75.              The number of columns of the matrix V.  M >= 0.
  76.  
  77.      V       (input/output) REAL array, dimension (LDV,M)
  78.              On entry, the matrix of right or left eigenvectors to be
  79.              transformed, as returned by STGEVC.  On exit, V is overwritten by
  80.              the transformed eigenvectors.
  81.  
  82.      LDV     (input) INTEGER
  83.              The leading dimension of the matrix V. LDV >= max(1,N).
  84.  
  85.      INFO    (output) INTEGER
  86.              = 0:  successful exit.
  87.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  88.  
  89. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  90.      See R.C. Ward, Balancing the generalized eigenvalue problem,
  91.                     SIAM J. Sci. Stat. Comp. 2 (1981), 141-152.
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.